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OO ' Abstract. 

3^ ' In the context of Brans-Dicke scalar tensor theory of gravitation, the cosmological 

^vq , Friedmann equation which relates the expansion rate H of the universe to the 

various fractions of energy density is analyzed rigorously. It is shown that Brans- 
5^. Dicke scalar tensor theory of gravitation brings a negligible correction to the matter 

C^ I density component of Friedmann equation. Besides, in addition to 17^ and Q.m in 

standard Einstein cosmology, another density parameter, fJ^ , is expected by the theory 
inevitably. Some cosmological consequences of such non-familiar case are examined as 
.,j. far as recent observational results are concerned. Theory implies that if Q.^ is found 

^ i to be nonzero, data can favor this model and hence this theory turns out to be the 

OQ \ most powerful candidate in place of the standard Einstein cosmological model with 

■^y ' cosmological constant. Such replacement will enable more accurate predictions for the 

^b , rate of change of Newtonian gravitational constant in the future. 
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Recent observational data have strongly confirmed that we live in an accelerating 
r^ \ universe [1] and have made it possible to determine the composition of the universe 

[2]-[l]. According to these observations, nearly seventy percent of the energy density 
in the universe is unclustered (dark energy) and has negative pressure by which it is 
driving an accelerated expansion [5]- [8]. Furthermore, the energy density of the vacuum 
is much smaller than the estimated values so far. By itself, acceleration seems to be 
much more understandable in the context of general relativity (cosmological constant) 
[9] and quantum field theory (quantum zero point energy); however, the very small 
and non-zero energy scale implied by the observations is not quite comprehensible. 
Because of these conceptual problems associated with the cosmological constant [TOj - 
|13j . alternative treatments to the problem have been produced and they are being 
used widely in the literature nowadays [14]-[l7j. For more detailed explanation about a 
number of approaches proposed so far and recent progress made towards understanding 
the nature of this dark energy see [18]. In some of these treatments, mostly, a scalar 
field is considered together with a suitably chosen V{(j)) to make the vacuum energy 
vary with time. 

To get a model in which the current value of the cosmological constant can 
be expressed in a more natural way; namely, without need of any fine tuning, in 
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the literature, there exist a number of studies on accelerated models in Brans Dicke 
theory [19]-[25]. For example, Sen et al [26] have found the potential relevant to 
power law expansion in Brans-Dicke (BD) cosmology whereas Arik and Qalik p^7J have 
shown that BD theory of gravity with the standard mass term potential (l/2)m^0^ 
is a beneficial theory in both explaining the rapid primordial and slow late-time 
inflation. In this regard, we have chosen the underlying theory as a scalar tensor theory, 
especially, BD scalar tensor theory of gravitation since scalar-tensor theories are the 
most serious alternative to standard general relativity. The theory is parameterized by 
a dimensionless constant cj, as c<j — i> oo BD theory approaches the Einstein theory j28j . 
In BD model, Lagrangian is defined by as in the following form 
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where the dimensionless uj is the only parameter of the theory, and Lm is the matter 
Lagrangian. According to our metric convention, (H ), Equation ([1]) turns out to be 
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In particular, it is expected that (p (t, x) is spatially uniform and evolves slowly only 
with cosmic time t such that ip (t, x) —>■ ^{t). As another point, the second term on 
the right hand side of ([1]) appears to be kinetic term of the scalar field but it is in 
an unlikely form, since the presence of the (p~^ which seems to indicate a singularity, 
and the presence of the coupling constant in multiplicative form is undesired. However, 
whole term can be transformed into the standard canonical form by re-defining the 
scalar field ip by introducing a new field 0, and a new constant e in such a way that 

V = ^.^' (3) 
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where e 

In this new form, BD Lagrangian is redefined as 

1 ,o„ 1 
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where the signs of the non-minimal coupling term and the kinetic energy term are 

properly adopted to (H ) metric signature in such a way that as g^^ ~ rf^, the 

kinetic term, \g^'" d^cpdycj) becomes |0^. Since the definition of a; is not changed, the 
limit uj approaching infinity is the same in both cases. In this limit ip remains constant 
whereas, as one can see from the above relation the limit of 0^ is singular. Since the 
matter part of lagrangian density Lm does not contain neither nor uj it does not 
take any part in this transformation and therefore is unchanged. Present limits of the 
constant u based on time-delay experiments [29]- [31] require uj > 10^ ^ L Besides, 
since these theories have been found as the low energy limit of string theory and they 
provide appealing models for inflation, scalar tensor theories enable an interesting arena 
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where the standard model can be tested. Hence, in this work, we aim to calculate the 
corrections, in the context of BD cosmology, to the famous Friedmann equation 

which relates the expansion rate H = a/a of the universe to the energy density. The 
fractional density parameter, Qi, is defined as the ratio of the energy density to the 
critical energy density which is a special required density in order to make the geometry 
of the universe flat. The Friedmann equation is used for fitting the Hubble parameter, 
H, to the measured density parameters {Q\, Qr, Qm) of the universe in such a way that 
Qa+ flR+ ^M = 1- According to recent observational results for the present universe, 
we have ^a — 0.75, Qm — 0.25, Qr ^ [32]. In the light of these values, one can 
conclude that the universe is mostly filled with non-baryonic matter and it seems that 
this non baryonic matter is responsible for the expansion of the universe solely. In the 
context of (BD) theory [33] whose scalar potential term consists of only a mass term 
and matter fields, the action in the canonical form is given by 
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The signs of non-minimal coupling term and the kinetic energy term are properly 

adopted to (H ) metric signature. As in the cosmological approximation, is 

spatially uniform, but varies slowly with time. As long as the dynamical scalar field (^ 
varies slowly, Gg//, the effective gravitational constant, is defined as G~}f = ^0^ by 
replacing the non-minimal coupling term g^0^ R with the Einstein-Hilbert term jq^q—R 
where R is the Ricci scalar in Einstein relativity. In natural units where c = h = 1, we 
define the Planck-length, Lp = uj/27r(f)Q where 0o is the present value of the scalar field 
(p. Thus, the dimension of the scalar field is L~^ so that the dimension of Geff is L^. 
Lm is decoupled from (p as was assumed in the original BD theory. Hence, considering 
does not couple to Lm as a matter field, we may consider a classical perfect fluid with 
the energy-momentum tensor T^ = diag (p, —p, —p, —p) where p is the pressure. The 
gravitational field equations derived from the variation of the action ([6]) with respect to 
the Robertson- Walker metric are 
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where k is the curvature parameter with k = —1, 0, 1 corresponding to open, fiat, closed 
universes respectively and a{t) is the scale factor of the universe (dot denotes d/dt), 
M denotes everything except the field. The right hand side of the equation ([9]) is 
set to be zero due to the assumption that the matter Lagrangian Lm is independent of 
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the scalar field 0. Since the matter lagrangian affects the first two equations ([31HD and 
when R is solved from the first two equations and substituted into the (j) equation ([9]) 
this equation depends on the trace of the energy momentum tensor T^^. Indeed this is 
how the BD field equations are sometimes written: 
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But in any case the system of equations and the solutions to the system of equations 
are independent of this phenomenon. Instead of working with the field equations (HEl) 
stated in terms of 0(t), a{t) and their derivatives with respect to the cosmological time 
t. We define the fractional rate of change of as F (a) = 0/0 and the Hubble parameter 
as H (a) = d/a, and rewrite the left hand-side of the field equations ([7119]) in terms of 
H{a), F{a) and their derivatives with respect to the scale size of the universe a (prime 
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From these three equations it can be shown that the continuity equation for the matter- 
energy excluding the BD scalar field is also satisfied with the help of the equation 

(in 
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and hence, instead of considering the p equation (1121) solely as one of the dynamical 
equations to be satisfied, we choose continuity equation in addition to the density 
equation and the equation to be satisfied in any cosmological case we want to explain. 
That is because once the continuity equation is satisfied than p equation must already 
be satisfied automatically provided that a is nonzero. To eliminate the dependence 
in (1 11 p . we take the time derivative of both sides of the p equation and after some 
rearrangements, we get flTT]) purely in terms of H{a), F{a), p{a) and their derivatives 
with respect to a. 
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After rewriting flT3l) in the following form 
ZaHH' - 2ujHaF' = -6H'^ 



3k 
2ujF'^ + QuHF - — + 2ujm 



(15) 
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we solve (1151 [ISD for H', F' and get the general form of the solution in the sense that 
once the curvature constant k and energy density in terms of a is given than H and F 
can be solved from the following equations: 
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Hence, in the present epoch, to discover how the Hubble parameter H changes with the 
scale size of the universe a, we assume that the present universe is mostly fiat and it 
necessarily obeys the pm = equation of state . Using flT^ . we find that the energy 
density p evolves with a in the same manner as in standard Einstein cosmology when 
the universe is solely governed by matter. 
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(19) 



where C is an integration constant. Setting k = and inserting this energy density into 
dm [IE]), we get the following form of the equations to be solved: 

-1 
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their derivatives with respect to u 

dH _ 1 

du aH{2uj + 3)u 
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we rewrite ( !20l [2T!) in terms of H{u), F{u) and 
[3 (2 + uj)H'^ + 2uj {uj + 1)F^- 2ujHF - 2uj{uj + l)m'^] (22) 



-H^ - (5c; + 6) F^ - 6(1 + uj)HF + um^ 



(23) 



Since these coupled equations are hard enough to be solved analytically for H and F, 
our approach is to determine a perturbative solution in which both H and F vary about 
some constants H^o and F^o respectively: 

H = H^ + Hiu + H2U^ + ... (24) 

F = Foo + FiM + F2u2 + ... (25) 

where Hoo, F^o, Hi, Fi, a, are all constants to be determined from the theory and from 
fitting the Hubble parameter, H, to the measured density parameters (^a, Qr, Qm) of 
the universe via Friedmann equation. Plugging this perturbative solution into 
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and keeping only the zeroth, first, second order terms of u and neglecting higher terms, 
we end up with two sets of solutions in the zeroth order 

v/(6cu2 + 17a; + 12)' '"2(^ + 1) ^ ' 

and 

H^ = ^^^^-F^ = Ih^. (27) 

Comparing the first order terms of m, on the other hand, provides two linearly dependent 
equations for which the only possible solution is the trivial solution of iifi = and Fi = 0, 

{[6{uj + 2) - a{2uj + 3)] H^o - 2ujF^} Hi + [-2ujH^ + 4uj{uj + 1)F^] Fi = (28) 

[-3/foo + Q{uJ + l)Foo] Hi + {[6 {io + I) - a{2uj + 3)] H^ + 2(5cu + Q)F^} Fi = 0. (29) 

Since the solution in which Hi and Fi are nonzero is much more plausible for our aim, 
the coefficient matrix is properly constructed from (128^ 1291) and its determinant is set 
to be zero to get the value of a for which Hi and Fi need not be zero simultaneously. 
We get two different a values 

a = 3 + 30 

l + UJ ^ ' 

and 

a ~ \fuj (31) 

corresponding to the solution sets ( l26l) and ( l27j) respectively. In this regard, we note 
two things here: 

• Concerning the solution of H we seek for, the solution (130|) is much more precious 
than the solution ( 1311) which approaches to infinity as lo becomes infinitely large. 
On the other hand, in the same limit, (|30|) gives a = 3 which is the well known 
term in a matter dominated universe solution of standard Einstein cosmology. 

• The correction factor 1/(1 + 0;) in the solution (l30l) is solely coming from the exact 
solutions of the field equations of BD theory. 

Two linearly dependent equations are available when one compares the second order 
terms of u\ 

{[6(cj + 2) - 2a(2cj + 3)] i^oo - 2ujF^} H^ + [4cj(cj + \)F^ - 2uH^\ F2 
= [a(2cj + 3) - 3(cj + 2)] Hi - 2tu{uj + l)Fl + 2ujHiFi (32) 

[3i^oo - 6(cu + 1)FJ\ H2 + {[2a{2uj + 3) - 6 (cu + 1)] H^ - 2{huo + Q)F^} F^ 

= - ^Hl + Fl + [6 (cu + 1) - a{2uo + 3)] HiFi. (33) 

Letting a = 3 + l/{uj + 1) and F^ = H^/2{lo + 1) in ([321 [33]) gives H2 and F2 only in 
terms of if 00, -^i, -Fi; 

^ ^ [-(3c^^ + 8cu + 6)^1^ + 2u{u + ifFl - 2o;(a; + 1)^1^1] 
'~ (3cu + 4)(2cj + 3)i7oo 
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{3u + 4){2uj + 3)Hoo ■ ^^^^ 

Hence, with these perturbation constants found from theory, we can express H and F 



as 



ff = i/. + i/.i^)''^ + i/.©' 
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where ( l26ll gives 



H^ = [2{uj+ 1) v^m] /v/(6a;2 + 17c; + 12) 



(36) 
(37) 

(38) 



and 



Foo = (v^m) /y/{6uj^ + nuj + 12). (39) 

To proceed one step further, we rewrite the standard Friedmann equation ([5]) with the 
extra term having new density parameter, Q^, 
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and we fit all theory parameters to the observational density parameters; 
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Hl = Hl + 2H^{H^ + H^) + Hi 
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(41) 
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With these relations above and the constraint Vt^ + VLm + fi^ = 1, we can express 
theoretical parameters Hi., H2 in terms of the observational density parameters f^A, VLm 



and Hr^ 
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Now at this stage, investigation of two cases for ^a can be meaningful: 



(45) 
(46) 
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We first set ^a — 0, ie; Qm + ^a = 1 consistent with today's universe density 

compositions and H2 7^ 0. By using recent observational results on density parameters 

^M ^ 0.27, ^A ^ 0.73 [35]- [36] together with (iS]) and ([461), we determine 

23 
H, = ^^T^^oo ^ 0.15/7oo, (47) 

H2 = -^ - _(H^ ^ _o.om^ < 0. (48) 



'■00 



Now, if (1361) is satisfied for H = Hq, we get 

H^ ~ 0.88iyo (49) 

where Hq is the present value of the Hubble parameter [31] and we may estimate 
mass of BD scalar field, m, for a fixed Hq as m '^ 10~^Hq by using the relation 
Hoc ^ 0.88i7o ^ 0.82m cji/2 for cj ^ 00. Furthermore, by using §^, ([35]), dSD, (jIH]) 
simultaneously, we get Fi ^ 0.08i7oo /^, -^2 — — 0.04ifoo /i^ for cu -^ 00. Remembering 
that Foo — Hoo /2ti;, we may say that F^o in (157|1 is the dominating term in today's 
universe. This shows us that similar to the expansion rate of the universe H, the rate of 
change of the Newtonian gravitational constant has approached the asymptotic regime. 
On the contrary, when we set ^a 7^ but ^a > ie; Qm + ^^a + ^a = 1 together 
with H2 — 0, from (H6l) . we immediately get the relation 

fiM = 2V^(l-V^), (50) 

and via this relation, we are able to estimate Qa and Q as in the following portions if 
the recent measured observational result on matter density parameter Qm — 0.27 [35] 
is kept fixed, 

^A ~ 0.71 

n^ ~ 0.02. (51) 

To figure out the reliance of such theory based density parameters and to show how this 
model can pass the observational constraint of the SNIa data [37] , we have compared our 
model of Hm = 2v^nA (l — a/Ha) proposing ^2 ~ 0.02, Qm — 0.27 and ^a — 0.71 with 
LCDM model of ^a = 1- ^m proposing Qm = 0.276 ± 0.026 and ^a = 0.724 ± 0.026 
according to wmap+SNIa data (ref). From such comparison, it can be seen that in an 
asymptotic regime, although ^a 7^ 0, SNIa data do not exclude our model as far as the 
range of uncertainties proposed by WMAP data team are valid. Besides, under this non- 
familiar case, we get other theory parameters Fi ^ 0.2Hoo /V^, -^2 — —7 x 10^^-ffoo /^ 
as cj ^ 00 . Remembering that F^o — H^o /2ti;, we see that Fi is the dominating term 
in (!37ll . namely, the rate of change of the Newtonian gravitational constant (Gn/Gn) 
has not approached to the asymptotic regime yet. However, theory predicts that when 
the size of the universe exceeds a ^ 0.6u;^/^ao, then the term F^ will become dominant 
so that asymptotic regime will be satisfied for f Gjm/Gn ) . 

As a last case, we investigate the possibility that Q^ < theoretically, ie; Q\ + 
Qm > 1- According to SNIa data, this case is important in the sense that if theory 
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assures that Q^ < then this will imply that data favors this model instead of standard 
Einstein cosmo logical model with cosmological constant. However, when we use f p^ and 
( 143|) simultaneously to make Q^ < we see that for Fi = Hi/2{uj + 1), H2 attains its 
most negative value of H2 = —H^/2H^ and for this value of H2-, ^^ exactly is being 
equal to zero instead of being equal to negative value. With this result we note that 
BD theory can not be forced to have Vt^ < so that data can favor this model instead 
of standard Einstein cosmological model with cosmological constant. 

Hence, in the light of these examinations on i7^, we may conclude that measurement 
of VL^ will be important in two respects. Firstly, if Vt^ is found to have positive value, this 
will indicate that data can favor this model instead of standard Einstein cosmological 
model with cosmological constant. Secondly, if that is so, making much more accurate 
predictions for the rate of change of Gn could be plausible. 

We would like to thank Burak Kaynak for his computational aid which provided us 
to be faster in seing the results of the work and Prof. Engin Arik for her valuable 
suggestions and contributions to the paper. This work is supported by Bogazici 
University Research Fund, Project no:06B301. 
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